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Correlated systems with hexagonal layered structures have come to fore with renewed interest 
in Cobaltates, transition-metal dichalcogenides and Gdl2- While superconductivity, unusual metal 
and possible exotic states (prevented from long range order by strong local fluctuations) appear 
to come from frustration and correlation working in tandem in such systems, they freeze at lower 
temperature to crystalline states. The underlying effective Hamiltonian in some of these systems 
is believed to be the Fahcov-Kimball model and therefore, a thorough study of the ground state of 
this model and its extended version on a non-bipartite lattice is important. Using a Monte Carlo 
search algorithm, we identify a large number of different possible ground states with charge order 
as well as valence and metal-insulator transitions. Such competing states, close in energy, give rise 
to the complex charge order and other broken symmetry structures as well as phase segregations 
observed in the ground state of these systems. 

PACS numbers: 71.45.Lr, 64.75.Gh, 71.30.-fh 



I. INTRODUCTION 

Systems like transition-metal dichalcogenidesi"— , 
cobaltatesi, Gdl-^^ and its doped variant Gdl2HjL 
have attracted considerable attention recently as they 
exhibit a number of remarkable cooperative phenomena, 
like valence and metal insulator transition, charge 
and magnetic order, excitonic instability^ and possible 
non-Fermi liquid states^iS,. The underlying geometric 
frustration of the triangular lattice, coupled with strong 
quantum fluctuations, gives rise to a large degeneracy 
at low temperatures and competing ground states close 
by in energy. A consequence of this is a fairly complex 
ground state magnetic phase diagram— and the presence 
of soft local modes strongly coupled to the conduction 
electrons''. 

These systems are characterized by the presence of 
localized and itinerant electrons confined to the two- 
dimensional triangular lattice. The electronic coherence 
along the perpendicular direction is negligible owing to 
the presence of large intervening ions (like Iodine in 
Gdl2) between the layers of rare earth or transition metal 
ions. It has been suggested recently that these corre- 
lated systems may very well be described by an effec- 
tive Falicov-Kimball model (FKM) on a triangular lat- 
tice. The nearly degenerate manifold of d-states close to 
the Fermi level, seen in LSDA calculations'", breaks down 
into a narrow, core-like state pushed below the Fermi 
level and an extended state crossing it, due to the strong 
local correlations in the d-manifold^ . The situation is 
similar to the well-known FKM^, if electrons in this nar- 
row, localized band are assumed to posses a very large 
effective mass. In the usual FKM one considers only 
the kinetic energy of the itinerant electrons and the local 



Coulomb interaction between the itinerant and localized 
electrons. This model has a long history and has been 
used to study metal insulator transition, magnetism in 
correlated systems and in the version where the inter- 
action is attractive (the sign of the interaction is irrele- 
vant in a particle- hole symmetric case) , it can be used to 
model crystallization in solidsiSi. Most of the studies are 
mean-field in nature. There are only a few exact results 
which are almost entirely on bipartite lattices — and 
are in lower dimensions'^. 

A very important question, in this context, is the na- 
ture of the ground states in such a system. Most of the 
systems described above show charge and orbital order 
at low temperatures. The layered material Naa;Co02, 
which shows superconductivity on hydration, reveals a 
\/3 X V3 superlattice'* at low temperatures due to spa- 
cial modulation of charge/orbital densities. Transition 
metal dichalcogenides mostly freeze into charge density 
wave statesiS, both commensurate and incommensurate 
(on doping some of them become superconductors at low 
temperature). Superlattice spots have also been seen in 
Gdl2 in X-ray scattering at low temperatures. Modelling 
the low-energy physics of some of these materials in an 
effective FKM naturally leads to such ordered states at 
low temperature. The FKM is known to support var- 
ious kinds of ordered or phase segregated'^' ground 
states on a bipartite lattice. Some of these ordered struc- 
tures appear in the ground state of a triangular lattice 
as well'^. There are hardly any rigorous results for the 
ground state of FKM on a triangular lattice. Most of 
the results are perturbativ o'^'^' , valid for large U. The 
combined effects of correlation and geometric frustration 
on the ground state phase diagram have not been investi- 
gated extensively using numerical methods on finite size 
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lattices either. With the appearance of charge ordered 
states as the low-energy broken symmetry ground states 
in the systems described above, it is imperative to look 
for the solutions of the FKM at low temperatures. In this 
work, therefore, we study the FKM in its spinless version 
numerically for all ranges of interaction, and explore the 
ground state charge ordering on a triangular lattice. 

It may be noted that in some rare earth compounds 
(specially the mixed- valence compounds), the rare earth 
ions with two different ionic configurations /" and /"~^ 
have different ionic radii. The contraction of d-orbitals 
in ions having /"~^ configuration, for reduced screening 
of nuclear charge, leads to increased localization of them. 
This implies that the d-orbital overlap between nearest 
neighbors depend on local /-electron configurations of 
neighboring ions, resulting in a correlated hopping of d- 
electrons. Such a correlated hopping term also appears 
in the first principles calculationi^ of the tight binding 
Hamiltonian and is usually neglected in Hubbard type 
models for being smaller compared to the on-site term. 
Its significance in superconductivity has been emphasized 
alread y^^'^" and in the context of FKM, it has been con- 
sidered by several authors^lr— using various approxima- 
tions. We, therefore, extend the model to include a cor- 
related hopping term in the Hamiltonian for the itinerant 
electrons. 

H = -Y.^U, + ^i5,,)d\d, +EfY^ flf. 

Here d\,di{jl , Ji) are, respectively, the creation and 
annihilation operators for itinerant d-electrons (localized 
/-electrons) at the site i. The first term in Eq.(l) is 
the kinetic energy of d-electrons on a triagular lattice: 
only nearest-neighbor hopping is considered. The second 
term represents the dispersionless energy level Efoi the 
/-electrons while the third term is the on-site Coulomb 
repulsion between d- and /-electrons. The last term, a 
correlated hopping term, is an extension over the usual 
FKM as discussed above. In the context of Gdl2 (or the 
transition metal dichalcogenides) , the /-electrons repre- 
sent the localized (or valence) band below the Fermi level 
while d-electrons cross the Fermi level in the effective low- 
energy model alluded to earlieri^i^i. 

II. METHODOLOGY 

The Falicov-Kimbll Hamiltonian (1), conserves the lo- 
cal /-electron occupation number nf,i = fjfi owing to 
its local U{1) gauge invariance in the absence of / — d 
hybridization. Therefore, [hf^i,H] = and uji = fi is 
a good quantum number taking values only or 1. The 
gauge symmetry also implies^^ interband excitonic aver- 
ages of the type < d\fj > are identically zero at any finite 



temperature (implying an absence of homogeneous mixed 
valence) and the /-electron level remains dispersionless 

iififj) — 0' fo'^ * 7^ j)- ^ non-bipartite lattice as we 
consider here, there is no particle-hole symmetry either. 
The local conservation of /-electron number implies that 
the Hamiltonian may be written as, 

H = Y. (^)dt dj +EfYl^^ (2) 

where %(a;) = [-ty + t'^ji^^i + uJj)] + {UuJi ~ ^i)5i.j. 

We set the scale of energy with t<y > = 1. The Hamil- 
tonian Eq.(2) represents a non-interacting d-electron 
moving in an annealed disordered background of the /- 
electrons^S. The interactions among the electrons is, 
therefore, kinematical in nature, coming entirely from 
the Fermi statistics for a given /-electron configuration. 
It suffices therefore to obtain the spectrum of this Hamil- 
tonian for different configurations w = {^1,^2, . . • ,WAf} 
of /-electrons by numerically diagonalizing the Hamilto- 
nian over different uj and annealing over the configura- 
tions. We perform this on a triangular lattice of finite 
size with periodic boundary conditions (PBC). Although 
this entails a tremendous simplification over solving a 
full interacting quantum Hamiltonian, even for a rela- 
tively small lattice size (in this work we considered L x L 
lattice with L = 12 and checked the results at larger L 
values in test cases), the number of different /-electron 
configurations is exponentially large and it may not be 
feasible to explore all the configurations within a reason- 
able computer time. For example, even for N=36 and 
Nf — 18, the number of configurations Ncj^^ ^ 10^'^. 

Therefore, one needs efficient search algorithms and 
explore only a relatively small subset of the entire /- 
electron configurations. We have used a Monte Carlo 
sampling method to achieve this goal. We work at half- 
filhng, i.e., Nf + Nd = N where Nf, Nj. are the total 
number of /- and c?-electrons and N = L'^ is the number 
of sites. For a lattice of N sites, the basis is chosen as 
4 |0 >,•••, djv |0 > and H (in Eq.(2)) is now an iV x TV 
matrix for a fixed configuration w. The partition function 
is, therefore 

Z = [] ( E Tre-^"^i^^A (3) 

where the trace is taken over the d-electrons, and /3 = 
The trace can be calculated from the eigenvalues 
Xi (i = 1 ■ ■ ■ N) of the matrix h (first term in H (Eq(2)) 
quite easily, h is diagonalized by a unitary transforma- 
tion WhU = K, where X is a diagonal matrix with 
its diagonal elements. The unitary rotation gives the di- 
agonal basis (in which h is diagonal) represented by the 
eigenvectors w| |0 >,•••, wjy |0 >. Defining the operator 
rii — vjvi and writing n.i to denote the corresponding 
eigenvalues of rii, the trace above can be identified with 
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FIG. 1: n/ — t' phase diagram for ?7 = 1, ?7 = 3 and U = 5. In the U = 5 phase diagram there is a stripe phase exactly at 
71/ = 1/2 from t' = — 1 to 0.8 (shown by a horizontal line from t' = — 1 to t' = 0.8) 



Tre~^^ = ^ <nx---nN\e~'^^\ni---nN > 

= ^ < ni---nAr|e~^Si^'"''|ni---nAr > 

(4) 

This reduction of the exponential to a c-number is the 
essential simplification that makes the problem amenable 
to a classical MC calculation. The partition function can, 
therefore, be recast in the form 

i ^i=0,l j = l 

= E e-^^f-')f[{l + e-^'n (5) 

i Wi— 0,1 J 

An average of any operator operating on the d-electrons 
is then calculated using the usual statistical mechani- 
cal procedure < A >= Tr{Ae~l^")/Tr{e~^"). Indeed, 



the intensive part of the calculation involves finding the 
eigenvalues of h. The summations over the static vari- 
ables {uii} arc evaluated using a classical Monte Carlo 
simulation. In order to get the ground state, a simulated 
annealing is employed ramping the temperature down 
from a high to a very low value. The process, there- 
fore, involves the following steps: (i) For a lattice of size 
L {N = L^), choose a particular value of Nf (0 < Nf < 
N, and Nd ^ N - Nj). (ii) For a fixed N and Nf, we 
choose a random configuration w = {lui,lu2, . . . ,uin} (hi) 
Choosing values for t',U and Ef, we find the eigenval- 
ues Aj of H{w) and the corresponding total free energy 
F{oj) = —kT log Z. The chemical potential /x is used to 
fix the (i-clcctron number Nd- The corresponding total 
energy is then E{oj) = limT->.o (iv) Generate a new 

random configuration oj' and calculate the new energy 

(v) Find AE = E{w) - E{w'), compare s = e'^^'^'^ 
with a random number r„ (0 < r„ < 1): if r„ < 
m?'n(l.s). accept the new configuration uj' , else it is re- 
jected. The steps (ii) to (v) are repeated until conver- 
gence is reached at a temperature T. The temperature 
is then ramped down slowly to a very low (compared. 



4 



(a) 



U= l,t'=-l,n„ 



1/4 



U= l,t'=0.2, n = 1/4 



o«o*o»o«o*o« 

oooooooooooo 
o«o»o»o»o«o» 
oooooooooooo 
o*o«o*o«o«o« 
oooooooooooo 
o*o«o*o*o*o* 
oooooooooooo 
o«o»o«o«o«o* 
oooooooooooo 
o*o*o*o*o*o* 
oooooooooooo 



•o»o«o*o«o»o 

oooooooooooo 
•o«o»o»o«o»o 
oooooooooooo 
•o*o«o«o«o«o 
oooooooooooo 
•o*o*o*o*o*o 
oooooooooooo 
•o«o«o«o*o«o 
oooooooooooo 
•o*o*o*o*o*o 
oooooooooooo 



U= l,t'= 1.0, n,-= 1/4 

oooooooooooo 
oooooooooooo 
oooooooooooo 
oooooooooooo 
oooooooooooo 
oooooooooooo 
oo»«oooooooo 
•••••••ooooo 

••••••••oooo 

o*«**«**oooo 
o*******oooo 
oo*»***ooooo 



(b) 



U= l,t' 



1.0, n= 1/3 



oooooooooooo 
•o«o*«o«o*o« 
o«o*o*o*o«o* 
oooooooooooo 
o*o*o«o*o*o* 
••oooooooo»o 
ooo«o«**o*oo 
••o»oooooo«o 
ooooo«o»o»oo 
••o«o«o*oo*o 
oo«o»ooooooo 
•o*oo*o*o*«o 



U= l,t'= 1.0, n^.= 1/3 

oooooooooooo 
oooooooooooo 
oooooooooooo 
oooooooooooo 
oooooooooooo 
oooooooooooo 
oooooooooooo 



••••oooo*««* 

••••oooooooo 



U = 5, t'= 1.0, n = 1/3 



oooooooooooo 
••••••••oooo 

••••••••••o^ 

••••••••••oo 

oo^^^^^^^ooo 

oooooooooooo 
oooooooooooo 
oooooooooooo 
oooooooooooo 
oooooooooooo 
oooooooooooo 



(c) 



U=l,t'=-1.0, n =1/2 



•ooao^a^ooo* 
o^o^o^oo^^o^ 
oo^^^^oo^o^^ 
••oaoo^^^ooo 
•ooo^o»oo^«o 
o^o^^oo^o^oo 
•••oo^^^oo^o 
•oo^o^oo^^^^ 
o«oo«^oo^o»o 
oo^o^o^o^o^^ 
•••••oo^^ooo 
•ooooaoao^ao 



U = 5, t'=-1.0, n = 1/2 



f 

••oo^^ooaaoo 
•oo^^oo^^oo^ 
oo^^oo^^oo^^ 
o^^oo^aoo^^o 
••oo**oo*^oo 
•oo^^oo^^oo^ 
oo^^oo^^oo^^ 
o^^oo^^oo^^o 
••oo^^oo^^oo 
•oo^^oo^^oo^ 
oo^^oo^^oo^^ 
o««oo«^oo**o 



U = 5,t'=1.0, n =1/2 



oooooooooooo 
oooooooo^^^o 



••••ooo^^^^^ 

oooooooooooo 
oooooooooooo 
oooooooooooo 
oooooooooooo 



(d) 



U= l,t'=-1.0, n, = 2/3 



o^^oo^^^o*^^ 
oo**^*^o*^«* 
••••oo^oo^oo 

••••oo^^^^^o 

•ooo^^^^oo^^ 
•oooa^^^^o*^ 
••••••oo^^^o 

o^^^^o^o^^^o 
o**o«oo*^^*^ 
•••oo^o^^oo^ 
••••o^^^^^o^ 

•ooo^^^oo^^^ 



U= l,t'=0.4, n,, = 2/3 

•o^^o^^o^^o^ 
o^*o«*o«*o«^ 
••o^^o^^o^^o 
•o^^o^^o^^o^ 
o^^o^^o^ao^^ 
••o^*o^*o**o 
•o^^o^^o^^o^ 
o^^o^^o^^o^^ 
••o^*o^*o^*o 
•o^^o^*o^*o^ 
o^^o^^o^^o^^ 



U= l,t'= 1.0, nj.= 2/3 

•••oooo*^**^ 

•ooooooo^^^^ 
•oooooooo^^^ 
••oooooooo^^ 
•••oooooooo^ 
••••oooooooo 
•••••oooooo^ 



FIG. 2: The ground-state configurations for various values of U, t' and n/. Filled circles correspond to sites occupied by 
/-electrons and open circles correspond to unoccupied sites. 



again, to typical AE) value and at each step of T, the 
same routine is performed. Repeating this over a few 
times usually lead to a unique low-energy state to the 
lowest temperatures searched. Although this is still a fi- 
nite temperature state, the nearest excitation AE was 
checked to be higher than KT indicating that this is 
the likely ground state. This procedure converges to the 
'ground state' faster (even for the largest lattice sizes 
we studied) in contrast with other methods^^ and makes 
it possible to test our results in larger lattices. There 
are, of course, the usual problems of simulated anneal- 
ing that one needs to take care. While ramping down 
in T, it often gets stuck in a local minimum even when 
the ramping is made ever slower at lower T. Different 
kinds of moves and at times simultaneous moves involv- 
ing several sites (and at longer ranges) were employed to 
drive the system out of the quench. The procedure has 
been successfully used in the search for complex ground 
state stuctures in manganites^''. It allows us to work in 
a much smaller sample-space of the whole configuration 
space, and any physical variable may be obtained by av- 



eraging over this sample space. Similar approach was 
employed in the study of the phase transitions in FKM 
on a bipartite lattice earlie r^^'^^ . 



III. RESULTS AND DISCUSSION 

Perturbative results for large U indicate that to or- 
der 1/U, the FKM can be mapped onto an Ising 
antiferromagnet (AFM) in a magnetic field, H^ff = 

J2<ij> iEF^i^J + hiP- + Ef)J2t + constant terms, where 
Si — 2uJi — 1, Si = —1, 1. The Ising AFM state on a 
triangular lattice is frustrated and leads to large degen- 
eracies at low temperature as discussed above. It turns 
out that this frustration is liftecti^ in the higher order 
perturbation in 1/U. It is therefore quite intriguing that 
one would expect the effects of frustration to bear on the 
ground states as U and the chemical potentials are var- 
ied. Different regions of the phase space are, therefore, 
controlled by different dominant effects. We studied the 
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FIG. 3: U — t' phase diagram at the /-electron fiUing Uf = 
2/3. Symbols showing the phases have the same meaning as 
in Fig. 1 



ground state configurations and the possible valence tran- 
sitions as a function of a range of values of U, and filling 
n/ = Nf/N {Nd is constrained to — Nf) using the 
method outlined in the previous section. We look at the 
effect of correlated hopping of d-electrons on these phases 
and obtain the ground state configuration for several val- 
ues of t' € [— 1, 1] at a fixed U. 

In Fig.l we present the phase diagram in n/ - i' plane 
at three representative U values {U =1,3 and 5). We ob- 
serve that the ground state configurations of /-electrons 
for a fixed filling are significantly affected by the corre- 
lated hopping t' and undergo several phase transitions at 
lower U values. At larger U, the ground states remain 
nearly independent of t' . Various ground state phases 
encountered in Fig.l are mainly; (a) Regular or quasi 
regular phase, found around nf = j (Fig. 2a columns 1, 
2), where the filled sites are arranged in an almost regu- 
lar pattern, (b) Bounded phase, where regions of empty 
sites are unconnected (and always surrounded by filled 
sites), observed around Uf — ^ (Fig. 2c and 2d, first col- 
umn). There is no apparent spatial order in this state 
and it is found to shrink as U increases, (c) Hexagonal 
or nearly Hexagonal phase (Fig. 2d second column) where 
/-electrons form hexagonal structures, are observed at 
all investigated U values. They appear at values of nj 
around | and |. (d) Segregated phase, where domains 
of sites occupied by /-electrons are segregated from the 
unoccupied sites (Fig. 2 column 3 shows four such config- 
urations). There can also be multiple such domains (not 
shown here), (e) Stripe phase: the filled sites form diag- 
onal stripes (Fig. 2c second column). This phase is ob- 
served at higher values of U (above U ~ 5) at = 1/2. 
It is found to be stable at finite t' ranging from —1 to 
0.8 beyond which there is a phase segregation. Compar- 
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FIG. 4: (color online) (a) Density of states for three different 
phases at n/ = 2/3 as is changed. Bounded: U = l,t' = —1 
(red), hexagonal: U = l,t' = 0.4 (green) and segregated: 
U — l,t' = 1 (blue), (b) Density of states for two different 
phases at n/ — 1/4. Regular: U = l,t' = —1 (black) and 
U = 5,t' = —1 (violet) and segregated: U — l,t' — 1 (blue). 
Ef is the corresponding Fermi level in each case. 



ing our results with earlier results at <' = limit^, we 
find that there can be a phase segregation now at a lower 
value of U (for example, J7 ~ 1) as the parameter t' is 
tuned. The t' term would, therefore, play an important 
role in driving phase segregation even in the weakly cor- 
related systems. This could be understood by looking at 
the position of the localized /-electron level {Ef) with 
respect to the d-electron band. At large U, the c?-band 
is strongly renormalized and the modulation of the ri- 
band due to the t' term becomes insignificant. It is also 
noticed that t' = 1 is a special point as at this point 
hopping takes place only between pairs of sites having 
identical /-electron occupation (near-neighbour configu- 
rations {0, 0} and {1, 1}, though, occupying different re- 
gions of the Brillouin zone) . This enhances the formation 
of inhomogeneous phases where a pair of /-electrons (or 
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FIG. 5: (color online) d-electron densities are shown on each site for U — 1, t' = —1.0, Uf = 1/4; U = 1, t' — 1.0, Uf — 
1/4; U = 5, t' = —1.0, Uf = 1/4 and U — 1, t' = 0.4, n/ — 2/3. The color coding and the radius of the circles indicate the 
d-electron density profile. The dark dot indicates cji = 1 at that site. 



a pair of empty sites) tend to occupy neighboring places 
on the lattice as in Fig. 2 column 3. If the coherence of 
the /-electrons is restored, this would lead to the forma- 
tion of interband excitons involving near-neighbor sites 
(at finite / — d mixing). In the context of Hubbard model, 
where d (/)-electrons represent up (down) spins, the cor- 
related hopping would, therefore, tend to enhance the 
extended s-wave superconducting order parameter (OP) 
fluctuations. 

Shown in Fig. 2 are some of the ground state patterns 
mentioned above and their dependence on the correlated 
hopping parameter t' as well as the Coulomb correlation 
strength U. Fig. 2(a) (top panel) reveals the variety of 
ground state patterns at n/ = 1/4 starting from regular 
structures to phase segregated regions on changing the 
value of t' from — lto-|-latC/ = 1. This is also ob- 
served at Uf = 1/3 where on varying either U or t' one 
could go from a quasi-regular pattern to phase segregated 
regions. At Uf — ^ (panel (c)) we have shown three pat- 
terns corresponding to three different sets of values for U 
and t' (i.e., bounded, stripe and phase segregated). On 
changing the value of U from 1 to 5, keeping t' fixed at 
— 1, a transition from a bounded phase to a stripe phase 
occurs whereas on changing the value of t' from -1 to 
+1 at a fixed U = 1, we arrive at a phase segregated 
state. And finally in panel (d) we present the phases ob- 
served at Uf — 2/3, ranging from bounded, hexagonal 
to phase segregated regions at three different t'. It is 



clear from the figure that a non-zero t' facilitates phase 
segregation in the ground state as it favors simultaneous 
occupation of /-electrons (or simultaneous occupation of 
vacancy, depending on the value of t') in the neighbor- 
ing sites. Phase segregations have been the key to the 
physics of many correlated systems hotly pursued in the 
last decade^^ and this correlated hopping term seems to 
open up another route to the phase segregation scenario 
in certain correlated systems. A ground state phase di- 
agram in the t' — U plane could be drawn now, based 
on the structures obtained in Fig. 2. We show this at 
Uf — I in Fig. 3. One can identify the three different 
phases, namely, hounded^ hexagonal and segregated. The 
reasonably large region of phase segregations at the upper 
part of the phase diagram owes its origin to t' primarily 
as discussed earlier. At large J7, of course, the hexagonal 
phase dominates. At this filling, no other phases were 
found in the range of parameters studied. 

In order to identify if there is a metal to insulator tran- 
sition in the above phase diagrams, we look at the d- 
electron spectrum (Fig. 4) and observe the gap"^^ at the 
Fermi level at various values of U and t' in the corre- 
sponding density of states (DOS). In Fig. 4(a) is shown 
the DOS for three distinct phases observed at n/ = 2/3 
keeping U fixed while changing t' . Fig. 4(b) shows the 
DOS for different values of J7, t' in two different phases 
at rt/ = J. We find that the hexagonal phase dX U — h 
is insulating with a large gap, the bounded phase at low 
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FIG. 6: Hf ~ Ef phase diagrams for U = 1 and different values of t' — —1, —0.6, 0, 0.6 and 1 are shown in (a) (e) respectively. 



U is also insulating but with a small gap while the phase 
segregated regions appear to have a very small (or vanish- 
ing) gap. Although the charge excitation gap continues 
to increase with U in the insulating phases, a 'spin' exci- 
tation involving exchange between vacant and occupied 
/-electron sites has a lower energy scale ~ ^ at large 
U. In this context, we note the evolution of the spec- 
trum and the corresponding DOS in Fig. 4(b) at a lower 
Hf — 1/4. It is clear that the DOS in the regular phases 
{U = l,t' = —1 and U = 5,t' = —1) have two-peak 
like structures above and below Fermi energy Ep (in the 
phases with a clear gap at the Fermi level) and the gap 
at the Fermi level increases with U. The band spreads 
out towards higher energy as expected. But the segre- 
gated phase {U = l,t' = 1) here is observed to have an 
almost gapless spectrum. Whether a very tiny gap still 
exists is not clear from this finite-size calculations, but 
within our numerical accuracy it appears to have closed. 
Such observations are reported earlier, albeit on a square 
lattice away from half-filling-'^'^ . It is possible that due 
to the segregation of /-electrons in one part of the lat- 
tice, d-electrons find a percolative path in the rest of 
the lattice and move freely from one end of the lattice 
to other without appreciable effects of correlation in the 
segregated phase. 

It is clear that the degeneracies as well as the sepa- 
ration between the energy levels increase as U rises and 
lead to incompressible states at specific values of the fill- 
ing (or, equivalently, the chemical potential) in the ther- 



modynamic limit. The FKM is known to show first order 
transitions and consequent phase separations in the bi- 
partite lattices and it appears similar physics holds good 
in large part of the phase diagram for a triangular lat- 
tice as well. In order to glean a physical picture of the 
local electron occupancies in different regimes, we have 
drawn the density of the d-electrons on the real lattice 
(Fig. 5) where the radius of the circles are proportional 
to d-electron density at that site (a color coding is also 
used) . A change of t' from 1 to — 1 leads to an increase in 
d-electron density on the unoccupied sites (sites without 
an /-electron) as discussed above. For t' = 1 there is a 
reasonable weight of the c?-electrons in the occupied sites 
at smaller U. As U increases the wave functions spread 
out from occupied sites. Thus the competition between 
localization and itinerancy is clearly visible from these 
figures. 

The FKM is known to have transitions 1- involving 
states with different Uf values. We have looked at these 
transitions as a function of Nf (keeping Nf + Nd = N) 
depending on the relative position of the /-electron level 
{Ef) and the Fermi level of the d-electrons. If i?/ lies 
above the Fermi energy {Ep), then Nd — N and the /- 
states are empty. In the opposite limit, we are in a clas- 
sical mixed- valent regime. For a given set oi t', U and 
Ef, one can easily find out the values of Nf for which 
the ground state has minimum energy. We have plotted 
Uf in Fig. 6 as a function of Ef for a set of values of t' 
and U. 
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Figs. 6(a). ..6(e) show Uf — Ef phase diagrams for [/ = 1 
and correlated hopping t' — —1.0, —0.6, 0, 0.6, and 1.0 re- 
spectively. We observe that as t' increases from —1.0 to 
+0.6, the valence transition (i.e. n^-transition) (i) oc- 
curs at smaller values of Ef and (ii) becomes sharper. 
However, as shown in Fig. 6(e), the ny-transition becomes 
smoother again for 1.0. We observe that the transi- 
tion width (the range of Ej over which nj goes from 
1.0 to 0.0), decreases as t' increases from —1.0 to 0.6. 
This could be explained from the fact that at t' = —1.0, 
the position of Ef, where n/ starts decreasing (from its 
maximum value 1.0) is located far below the centre of the 
d-band. As t' increases towards zero the position of Ef, 
at which Uf starts decreasing, moves towards the centre 
of c?-band. For t' > 0.6, the corresponding position of Ef 
again moves below the centre of d-band. So the effective 
width of d-band is maximum for t' = —1.0, is minimum 
around t' — 0.5 and increases again for t' > 0.5. In 
Fig. 7(a), 7(b) and 7(c) we show the same for U = 1,5, 10 
with t' — —0.6, and 0.4. We observe here that the 
transition width decreases as the value of U increases. 
The value of Ef where major change in n/ occurs hardly 
depends on the value of U. For all values of U, the rt/- 
transition is relatively smooth for t' around —0.6 and this 
transition becomes steep for t' around 0.4. These features 
are easily understood from an analysis of the spectrum 
of d-electrons. 

As an interesting aside, we note that although the ex- 
tended FKM studied here has U > 0, the corresponding 
ground state phase diagram for J7 < is obtainable from 
the following symmetry: if {uJi} is the ground state for a 
particular /i, Ef and U, then {(D^} is the ground state for 
/i, —Ef, —t' and —U, where Ui = 1 — Wi. The negative U 
FKM was introduced by Lieb^^ as a model for the study 
of the formation of crystalline solids. In that context, the 
heavy /-electrons can be thought to represent the ions in 
a solid, and their ordered structure implies a crystalline 
arrangement of the lattice. 

It is in order now to discuss the numerical results ob- 
tained above in the context of the systems where such 
models are expected to be useful under simplifying as- 
sumptions, in certain regions of their phase diagrams. 
Fig.l gives a good account of the complexity of phases 
obtainable from an FKM on a triangular lattice. The 
model produces charge ordered states of varying period- 
icities and phase segregated states as filling, U and t' 
are changed, also observed in dichalcogenides, cobaltates 
and Gdl2. The hexagonal structure close to n/ = 2/3 is 
reminiscent^ of the superlattice structure in Na:rCo02 
(at X = 2/3 in Fig.3(c) of Ref. 4). The CDW in 2H- 
TaS2 is incommcnsurate^i, while it is commensurate in 
TiSe?^, though the superlattice periodicity is dictated 
by special relations between the nearly flat valence band 
and the conduction band^!^, rather than the geometry of 
the underlying lattice, since the CDW in TiSe2 is likely 
to be excitonic in origir^. Stripes in correlated systems 
tend to appear at special fillinga^l due to commensura- 
bility effects in bipartite lattices, while we observe its 



formation here (e.g., at Uf = i) even if the underlying 
lattice is non-bipartite. Stripes have not yet been seen 
experimentally in the systems we discussed, but appear 
to be a distinct possibility if carefully investigated. The 
phases revealed here clearly shows the intricate nature 
of the competing states alluded to earlier - correlation 
and frustration working together on a triangular lattice 
in unraveling such a rich phase diagram. This under- 
lines the need for a thorough experimental investigation 
of these systems to delineate both the ordered and inho- 
mogeneous phases as filling and parameters change. 

As discussed earlier, a common feature among systems 
like dichalcogenides —1^ and Naa;Co02 — is the competi- 
tion between a charge density order and superconductiv- 
ity, where slight changes in symmetry or a small pertur- 
bation may lead to a transition from one to the other. 
Off-diagonal long range order (ODLRO) is absent in the 
FK model unless of course the /-electrons become disper- 
sive. In the limit of infinitely heavy down-spin electrons 
(iyU — 0), the Hubbard model reduces to the (sym- 
metric) FKM. With any symmetry-breaking perturba- 
tion (pressure^, doping/hydration^r^ etc.), as i^j-.j, shifts 
away from zero, the charge order or phase segregation 
may disappear— and superconductivity could appear: s- 
wave for the attractive-C/ model and extended s-wave in 
the repulsive-C/ case, due explicitly to the t' terrain. The 
nature of superconductivity in both these systems should, 
therefore, shed considerable light on the underlying mi- 
croscopic mechanism operative here. Although there is 
a suggestion of high angular momentum superconduct- 
ing OP in doped dichalcogenides-, the symmetries of OP 
in hydrated cobaltate and doped dichalcogenides are not 
resolved yet and should be probed thoroughly. 

In Gdl2 a phase segregation appears to preempt the 
formation of a long range order. While phase segrega- 
tion has been shown to appear in FKM on bipartite lat- 
ticeaii'^, it requires a fairly strong U . Gdl2 falls in the 
intermediate coupling range and the phase segregation 
seen in Gdl2 at low temperatures even for a low U could 
be facilitated by the presence of t' as shown above. The 
correlated hopping term helps in tuning the balance be- 
tween states close in energy, and therefore, plays a signif- 
icant role in non-bipartite lattices where nearly degener- 
ate states proliferate at low temperatures. We note that a 
quantitative modeling in detail, including first principles 
band structure and Fermi surface topologies, is beyond 
the present model. Such a study is indeed very useful and 
left for the future. We believe our present study already 
sheds light on certain important aspects of the ground 
state phenomenology and motivates new experiments. 

In conclusion, we have studied the Falicov-Kimball 
model on a non-bipartite lattice and found several or- 
dered ground states. Extending the model to include a 
correlated hopping term leads to new effects: in particu- 
lar, it strongly facilitates phase segregation in the ground 
states even in the weak correlation limit. We also identify 
several valence transitions as i' and i7 are changed. These 
observations are quite relevant for the ground state or- 
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FIG. 7: Uf — Ef phase diagram for different values of correlated hopping t' = —0.6,0 and 0.4, at (a) 17 = 1 , (b) (7 = 5 and 
(c) U = 10. 



der and phase segregation observed recently in correlated 
systems with triangular lattices such as Gdl2, dichalco- 
genides and cobaltates. 
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